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Mechanical dissipation poses an ubiquitous challenge to the performance of nanomechanical de- 
vices. Here we analyze the support-induced dissipation of high-stress nanomechanical resonators. 
We develop a model for this loss mechanism and test it on SiaN4 membranes with circular and square 
geometries. The measured Q- values of different harmonics present a non-monotonic behavior which 
is successfully explained. For azimuthal harmonics of the circular geometry we predict that destruc- 
tive interference of the radiated waves leads to an exponential suppression of the clamping loss in 
the harmonic index. Our model can also be applied to graphene drums under high tension. 

PACS numbers: 85.85. +j, 42.50.Wk, 63.22.-m 



Nanomechanical resonators offer great potential for 
practical device applications that exploit their ultra-low 
mass and high frequencies^. Examples range from scaling 
scanning-probe force microscopy and mass-sensing down 
to the molecular scale to providing alternatives for radio 
frequency devices. These applications share the desirabil- 
ity of high mechanical quality-factor (Q) that, by virtue 
of narrow bandwidth, amounts to a better defined fre- 
quency thus enhancing performance. In turn, measure- 
ments of mechanical displacements with an imprecision 
below the standard quantum limit and the preparation 
of ultracold motional states have already been imple- 
mented with electromechanicaPSl and optomechanical 
systems^'. These breakthroughs foreshadow the possi- 
bility of realizing a "quantum optics" analogue involving 
a macroscopic mechanical degree of freedom which would 
set a new st age fo r fundamental tests and potential quan- 
tum devices^EEl. Once more mechanical dissipation, as 
determined by the Q-value, plays a critical role in such 
endeavors. 

Though the mechanical Q-value may in general be in- 
fluenced by various mechanisms, in a small suspended 
structure that is sufficiently clean and cold, internal losses 
induced by two-level nuctuator^HH anc j radiation of elas- 
tic waves into the substrate are likely to play leading 
ro i e ^i3iti51 i n this respect, the reduction of the design- 
limited "clamping loss" induced by the coupling to the 
substrate will allow the use of nanomechanical devices to 
probe the internal losses and quantify the fundamental 
contributions of the constituent materials. Furthermore, 
with the advent of the use of stressed silicon nitride mem- 
branes, the high Q-values of these devices have demon- 
strated that it is indeed possi ble to attain low mechanical 
losses in nanoresonatoraUBZOES 

In this letter, we present and test a model that cap- 
tures the energy loss that occurs due to elastic wave ra- 
diation at the periphery of these high-stress resonators. 
We show that this mechanism is significant in state of the 
art structures and is strongly influenced by interference 
effects. We compare the results of our model to mea- 
surements of the resonant modes of two configurations, a 
single "drum resonator" and a composite array of drum 



resonators that effectively realizes a square membrane 
[cf. Fig. 1 (a) and (b)]. We examine the harmonics of 
these structures and accurately account for much of the 
variation in the corresponding Q-values. Our analysis re- 
veals that certain types of modes are inherently resilient 
to clamping loss as a result of destructive interference of 
the radiated waves. Thus, we provide insight into res- 
onators that might be realized and yield better Q-values 
in the future. On general grounds, the fact that the rel- 
evant stress at the resonator-support contact scales at 
least linearly with frequency combined with the 3D na- 
ture of the support, lead to the naive expectation that 
the dissipation (1/Q) due to elastic-wave radiation should 
increase as one considers higher harmonics [cf. Eq. Q]. 
In dramatic contrast, we find that for the harmonics of 
a circular membrane the clamping loss is exponentially 
suppressed as the number of radial nodal lines increases. 

To derive an adequate model for the clamping losses, 
we adopt the phonon tunneling approach introduced in 
Ref. HH and start from the general weak coupling expres- 
sion for the dissipation 1/Q in terms of the "overlaps" 
between the resonator mode and the free modes of the 
substrate ( "support" /221 
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Here cr' R and u' R are the stress and displacement fields as- 
sociated with the normalized resonator mode, er^ and 
uf^ are the analogous fields for the continuum of sup- 
port free modes labeled by q [eigenfrequencies uj(q)], and 
p s and p R are, respectively, the densities of the sub- 
strate and resonator materials. In our setting the res- 
onator mode should satisfy clamped boundary conditions 
at the resonator-support contact area S while the un- 
perturbed support modes should satisfy free boundary 
conditions implying that only the second term in Eq. 
contributes. The substrate is modelled as a half-space 
that contacts the membrane resonator at its rim S - 
i.e. the underetched gap between the suspended structure 
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FIG. 1. (a) Micrograph of a single drum resonator (similar to the one used in our analysis) superposed with schematic diagrams 
of the different harmonics (n, m) depicting their nodal lines (the origin is set at the center of the membrane), (b) Idem for a 
square membrane resonator, (c) Dissipation 1/Q as a function of frequency for the different harmonics of a SiaN4 drum resonator 
(D = 14.5/jm, t = HOnm, and a = 0.90GPa). (d) Idem for a square membrane resonator (253.2/xmx253.2/imx0.0125/im 
with a = 0.87 GPa). Red plot: Measured values with an error of 10% for 1/Q — we ascribe the splitting of degeneracies 
observed for the square membrane [cf. (d)] to disorder. Blue plot: Least squares fit of our model to the measured 1/Q using as 
fit parameters an internal dissipation offset (1/Qint) and material properties of the substrate (E s — 148 GPa, p s = 3.75gcm -3 , 
and 1/Qi n t = 8.5 x 10 -7 for the square membrane, and E s = 323 GPa, 1/Qi nt = 4.6 x 10 5 for the drum). Green plot: 1/Q 
without the offset corresponding to the predicted clamping loss — the resulting limits for the Q- values of the drum and square 
membrane are shown, respectively, in (e) and (f). High-Q harmonics (n,l)\ n >o [(n, n)\ n >i] of the drum [square] are marked 
by circles; low-Q harmonics (n, l)| n >o of the square, by triangles. For the drum [square] all harmonics with frequencies below 
130 MHz [9 MHz] are included except (0,3), (3,2) [(7, 1)]. 



and the substrate is neglected when determining the sup- 
port free modea 2 ^. We assume the "high stress" regime 
t 2 /D 2 <C o/E R <C 1, where a is the tensile stress in the 
membrane, t its thickness, D its large dimension (diame- 
ter or side) and Er the Young modulus of the resonator 
material. This implies that bending effects arc negligible 
and one can use the classical wave equation adequate for 
a taut membrane^. Thus for the drum's eigenfrequen- 
cies we obtain [cf. Fig. 1 (a), (b)]: u> nm — 2C, nm cji/ D 
with n — 0, 1, . . . and m = 1,2,...; while the square's 
eigenfrequencies are given by: ui nm = n^/n 2 + m 2 Cji/D 
with n, m = 1, 2, . . . Here cr — \J a / ' pn is the phase ve- 
locity in the membrane, and £ nm is the mth zero of the 
Bessel function J n (x). In this context, the weak coupling 
approximation underpinning Eq. reads Lu nm t/cR <C 1. 

For the single drum we adopt support eigenmodes 



(r) (with i = Q,±l, ...) that have axial symme- 
try with respect to z [cf. Fig. 1 (a)]. These are related 
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to the plane wave eigenmodes u^(r) by 
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l{-i) n /V^]fl„dipe mv u^l / (f); where j = l,t,s labels 
the different types of relevant modes [i.e. longitudinal 
(I), transverse SV (t), and SAW (s) given that SH waves 
do not contribute] with velocities of propagation c 7 , and 



we use spherical coordinates for the incident wavevec- 
tor q(q, 9, ip) — q(sin9 cos (p, sin 9 sin <p, cos 9) [9 = ir/2 for 
7 = s and 9 < ir/2 otherwise]. We note that the stress 
<t' r corresponds to the variation with respect to equilib- 
rium and that the resonant wavevectors in the substrate 
satisfy w nm /c 7 t^Cl, so that we can neglect the variation 

of u^l j (f) across the thickness t. Thus, substitution of 
the support and resonator modes into Eq. ([!]) (cf. Ap- 
pendix IaI leads to 
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Here we introduce the dimcnsionless functions 
ui :1 ^ s {q,v s ) = 27r/ 7r/2 d6»sin6'|4 7;z (0,^)| 2 J ; 2 (qsin6l), 

ui rS (q,u s ) = 27r|z4°]. z (0, v s )\ 2 jf{q) and define 

?7 7 = Cfl/c 7 ~ \J crp s /E s pii — where the prefactors 
of order unity, which depend on 7, are functions 
of the Poisson ratio for the substrate v s . We note 
that |itg° 7 . z (0, v s )\ 2 solely depends on cos 9 and v s 



cf. Eqs. (A9)-(A12| 



In turn, for the square membrane an analogous proce- 
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dure detailed in Appendix [5] leads to: 

1 16irn 2 m 2 p R t ^ 3 _ / -^- 

Wnrn V™ + m z p s D , 

(3) 

with w"; m (g, t/ s ) = \u^l ;z (0,u s )\ 2 f nm i(q), w^ s (q,v s ) = 

J*' 2 d6sin6\uf^. z (0, v s )\ 2 f nm i(q sin 9), where I > and 
we introduce f nm i{x) = fi(wx)[Z nm i(x) + Z mn i(x)\ with 



Znml (*e) : 
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Here z„g(x) = 2n 2 — x 2 =F 2n\Jn 2 — x 2 and the func- 
tions /;(a;) are given by: fi(x) = [8 i - 2(-l)"J 4; (x) + 
{-l) l J 4l (V2x)}/(2 Sal x 41 ) for n + m even, and //(a;) = 
[6 0l ~ 2 sin 2 (/tt/2) J 2l (x) - cos (Ztt/2) J 2/ ( >/2x)] / (2 5w x 21 ) for 
n + m odd. Eq uation ([3] ) is only valid for the case 
mm{n,m} > rj s \/n 2 + m 2 which is satisfied for the res- 
onances studied here — note that material properties 
always imply r/i < i] t < rj s , and a <C E s implies tj s « 1. 

We proceed to compare the predictions of our model 
[Eqs. ([2|-([3|] with the dissipation measured in nanome- 
chanical membrane resonators (cf. Fig. 1). These res- 
onators are made of "stoichiometric" SisN4 deposited by 
low pressure chemical vapor deposition on SiOa 11 . The 
nitride has an inherent stress of 1.2 GPa, as measured by 
a wafer bow technique^, and a density pn — 2.7gcm -3 . 
After lithographic patterning to define access holes, the 
resonators are suspended by etching the underlying oxide 
through these holes, using buffered oxide etch (BOE) for 
the single drum and HF for the square membrane, and 
critical point dried. Thus a single access hole results in a 
circular drum geometry, while a square geometry is de- 
fined by a periodic square lattice of such holes (50 x 50 
separated by 5 /jm) . Given the small size of the holes 
(< 1 /im) compared with the typical mode wavelength, 
we neglect them in our model. For the square array, the 
same consideration applies to the hole separation so that 
we use a square membrane model with uniform thickness 
t = 12.5 nm given by the average over the array^ and 
side D = \[A — 253.2 pm, where A is the suspended 
area. For the single drum (diameter D = 14.5 pm) the 
use of a BOE etch implies that the thickness is uniform 
and equal to the nitride thickness (llOnm). 

The mechanical resonances of the structures are char- 
acterized under vacuum and room temperature condi- 
tions, using a technique described in Ref. 111! The res- 
onators are actuated using a piezo disc that vibrates 
the chip in the out-of-plane direction and the motion is 
detected via a 633nm continuous wave laser. Figure 1 
(c)-(f) compares the measured frequencies and Q-values 
of different harmonics for the two configurations, single 
drum and square array, with the predictions of our model. 
This comparison takes into account three issues: (i) the 



release of the resonator leads to a local deformation of 
the wafer that lowers the membrane's tensile stress with 
respect to the one in the nitride layer, (ii) in addition to 
clamping losses the resonator will be affected by internal 
dissipation, and (hi) the parameters for the half-space 
model of the substrate must be judiciously chosen. 

To deal with (i) we determine the membrane phase 
velocity cr from a suitable linear regression that uses 
as input the resonator size D, the measured frequencies, 
and their mode indices which can be identified from the 
frequency ratios between the harmonics and the funda- 
mental mode. We find an excellent correlation that yields 
cr = 576.8ms -1 (566.8ms -1 ) for the drum (square). 

Our model is in excellent agreement with the observed 
trends providing the internal dissipation channel (ii) is 
frequency independent, and can be just added as a fit 
parameter (1/Qint) to the calculated dissipation. To elu- 
cidate (iii) one needs to compare the wavelengths of the 
resonant "support" modes with the thickness of the Si 
wafer (0.5mm). For the square the resonant frequen- 
cies are in the MHz range resulting in wavelengths in Si 
(4-8mm) much larger than the wafer's thickness so that 
these modes are dominated by the properties of the un- 
derlying piezo and positioning system. Thus, we adopt 
v s = 1/3 and leave the density p s and Young modulus 
E s as fit parameters. On the other hand for the drum 
the resonances studied lie in the 100 MHz range so that 
the elastic wave radiation is determined mostly by the 
anisotropic properties of crystalline Si. For this case we 
adopt p s = 2.33 gem -3 and v s = 0.28, but leave E s as a 
fit parameter given the isotropic nature of our model. 

In both cases, drum and square geometry, we find a 
class of modes that consistently exhibit lower dissipation 
1/Q when compared to nearby modes [cf. Fig. 1 (c) and 
(d)]. Their measured Q remains approximately constant 
as the harmonic index is increased, leading to a growth in 
their fQ product that for the square reaches a maximum 
of 1.0 x 10 13 Hz for the (6,6) harmonic. These "special" 
classes of harmonics for the drum and square are, re- 
spectively, (n, l)| n >o and {n,n)\ n> i and correspond to 
the presence of nodal lines that intersect the periphery 
at evenly spaced points [cf. Fig. 1 (a) and (b)]. In con- 
trast, for the square geometry the modes (n, 1), (l,n), 
where two of the sides do not intersect any nodal lines, 
tend to exhibit smaller Q for comparable frequencies with 
fQ ~ 10 12 Hz. An intuitive heuristic understanding of 
these trends emerges from realizing that for low harmon- 
ics, with membrane wavevectors ~ n/D, the typical res- 
onant wavelengths in the substrate are much larger than 
D. Thus, for the special modes the clamping loss is sup- 
pressed [cf. Fig. 1 (e) and (f)] due to destructive interfer- 
ence between the waves radiated by the different equiva- 
lent segments of the periphery, defined by the nodal lines, 
which have alternating 7r-phases. Concomitantly, unlike 
the fundamental mode, these special modes are associ- 
ated to stress sources with vanishing total force. 

A quantitative grasp of these striking features can be 
gained by exploiting the smallness of the ry 7 underpin- 
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ning the aforementioned wavelength separation. For the 
drum, relevant harmonics satisfy the condition ?y 7 Cnm <C 
y/n+l which allows us to Taylor expand the Bessel func- 
tions in the w;. 7 yielding an approximation for Eq. ([2]) 
that implies the following^ 



Qoi 

Qnl 

Qoi 



2n 2 a R c 2 R uJoiUo(v s ) 

2n 



0.517- 




where <jr — put is the surface mass density of the 
membrane and Uq(v s ) = J2- y ( c t/ c -y) 3 'Uo yl (0,b' s ). Thus, 
the clamping-loss limited Q-values of modes (n, 1) ef- 
fectively grow exponentially — as the super-exponential 
factor plays a negligible role for relevant n, in sharp 
contrast to series of modes for which m is increased 
while n is kept constant. These exhibit a decrease of 
(Qciamp for increasing frequency^. On the other hand, 
for the square geometry analogous considerations imply 
for m ~ tl ^> £oi/27T?7 7 a rise in Q c iamp that is merely lin- 
ear, with the damping rate tending to a constant value, 
as the harmonic indices are increased with their ratio 
m/n fixed. In turn, for our setting given the magnitude 
of 1/Qi ri t all the high-Q modes present roughly constant 
Q- values. 

A comparison between the predictions [cf. Eqs. ^ and 
([3])] for the two geometries (with appropriate dimensions) 
also reveals that for "special" harmonics [(n, l)| ra >o and 
(n, n)\ n> i] with the same frequency and number of nodal 
lines the circular geometry always yields a higher Q. Fi- 
nally, one should note that the scalings, embodied in 
Eq. Q, for the Q-values in terms of p R /p s , E s /a, and 
D /t are completely general and independent of the shape 
of the boundary. These, directly imply that the /Qciamp 
product of a given harmonic is independent of D. Fur- 
thermore, typical parameters yield for the fundamental 
mode /Qdamp ~ 10 12 Hz, which is comparable to experi- 
mental values (cf. Fig. 1 (c) and (d), and Refs. I2"0l |2T|) . 

In conclusion, we find that the dissipation of different 
harmonics of a given membrane resonator exhibit a strik- 
ing non-monotonic behavior which can be understood in 
terms of how the mode-shapes of different harmonics in- 
fluence the clamping loss. We find classes of modes for 
which the measured Q remains approximately constant 
and substantially larger than for other modes with com- 
parable frequency, and explain this phenomenon in terms 
of destructive interference between the radiated waves 
leading to a strong suppression of the clamping loss. No- 
tably, our analysis implies that for modes (n, 1) of a cir- 
cular geometry, the damping rate due to elastic-wave 
radiation vanishes exponentially in n rendering them 
"asymptotically mute". Thus, for typical parameters, 
these azimuthal harmonics can be regarded as effectively 
clamping-loss free for moderate n (e.g. /Q c iam P ^ 10 17 Hz 
for n > 5 and thickness t < 100 nm). Our results 
are relevant to state-of-the-art dispersive optomechanical 



setup^212U anc i the model is also applicable to graphene 
nanodrums under tension^!. Finally, we highlight that 
the interference effects we have unveiled will also be op- 
erational for the flexural modes of rigid plates. 

IWR acknowledges financial support via the Nanosys- 
tems Initiative Munich. Work at Cornell was supported 
under NSF ECCS 1001742. 



Appendix A: Circular geometry 

In the case of the circular geometry, the axially sym- 
metric support eigenmodes used are related to the plane 



wave eigenmodes 1A ' (r) by: 
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Here 7 = I, t, s labels the different types of rele- 
vant modes [i.e. longitudinal (I), transverse SV (t), and 
SAW (s) given that SH waves do not contribute] and 
q(q,9,ip) = q(sin9 cos ip, sin 9 sin ip, cos 9) is the incident 
wavevector [0 = n/2 for 7 = s and 9 < tt/2 otherwise]. 
As the resonator material is prestressed, the stress cr' R in 
Eq. (1) is not the total stress but corresponds instead to 
the variation with respect to equilibrium. This, together 
with the validity of membrane theory)^ implies that to 
linear order in the displacement field we have 



dS ■ <r' 
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which we will use to simplify Eq. (1). Here we intro- 
duce the normalized resonator eigenmodes 4>r — ¥ 4> nm , 
that satisfy the 2D classical wave equation with <p nm = 
at the periphery of the membrane, and adopt cylindri- 
cal coordinates f=(r cos <j>, r sin <j>, z). These eigenmodes 
read 



0om (r,4>) = d~R 0m (r) , 

4>nm{r,4>) = \j^Rnrn(r) COS n(j) , 

(f>nm{r,4>) = J -Rnm(r) sin ncj) 



with 



Rnm (^") 



2J n (2( nm r/D) 



DJn+l (Cnm) 

and have eigenfrequencies U)r — > u> nm given by 
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Here n — 0, 1, . . ., m — 1, 2, . . ., and the harmonics with 
n > are doubly degenerate. 

As the resonant wavevectors in the substrate satisfy 



i/cyt <C 1 we can neglect the variation of 



(o) 



(f) 



across the thickness t (i.e. the z-dcpendence at S). This 
approximation and Eq. ( A2 ) imply that the support 



modes only enter into Eq. (1) through u'fl t (r)\ z =o- 
To determine the latter we exploit that reflection at the 
free surface preserves the tangential component of the 
wavevector implying 
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Then, from Eqs. ( Al ) and ( A7), using the Bessel integral 
Jn(x) = ± f^ e ^^-*smV7d<£, we obtain 
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We now deploy Eq. (|A2j) and substitute Eqs. (|A3|)-(|A5|) 



and (A8| into Eq. (1). Subsequently, we use that 
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where d 7 is the dimensionality (i.e. d 7 = 2 for 7 = s 
and d 7 = 3 for 7 7^ s), perform the substitution uj = c~ f q 

(for each 7), evaluate using J' n {Q nm ) = -J n +i{C,nm), 
and express cr in terms of cr and pr. Finally, integration 
over uj and <fi, summation over I (all terms vanish except 
l — ±n), and substitution of Eq. (A6) yields Eq. (2). 



The ^-displacements at the origin (0) only depend 
on cos 9 and v Sl and their absolute values squared are 
given by the following^ 
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for < u < Vl - a , (A10) 
for vT^a < u < 1 , (All) 
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Appendix B: Square geometry 

For the square geometry we exploit that the whole 
structure still presents reflection symmetries with respect 
to the x — z and y — z planes. These are associated, re- 
spectively, with the operators R x and R y acting on the 
space of solutions of the elastic wave equations. Thus we 



can use normal modes of the decoupled support (elastic 
half-space) and of the resonator (i.e. square membrane) 
that are eigenvectors of R x and R y . In the case of the 



support, for a given plane wave mode \u) 



,(0) 



(r) of the 



elastic half-space one can generate modes \u ++ ), ), 

I it |_), and \u ) with the desired reflection properties 

by the following symmctrization procedure 



= S, 



//// 



jj,R x \u) 



VRy\u) 



/M/R X Ry 
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with n, v = ±, which enforces 



(B2) 



Naturally a complete basis is obtained by taking q in the first quadrant, i.e. < tp < Tt/2. As for the drum: (i) 
the smallness of t compared with the resonant wavelengths in the substrate implies that we can neglect in Eq. (1) 
the z-dependence of the symmetrized support modes at S, and (ii) only their z-components are relevant. These 
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^-components read 



uf^ ++ . <z (x,y,0) = 2 4t; 


z (0) cos(qx sin 9 cos ip) cos(qy sin 9 sin tp) , 
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, (0) sin(<7a; sin 9 cos </?) cos(qy sin 6* sin . 
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In turn, for the resonator modes we obtain 
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for n, m even — ^ antisymmetric case , 


(B8) 


2 n7r m7r 


for n odd and m even — > mixed symmetry case , 
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for n even and m odd — > mixed symmetry case , 
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with eigenfrequencies 
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Naturally only modes with the same reflection symmetries are coupled via Eq. (1). We first treat the symmetric 



modes and then briefly outline the straightforward extension to the other cases. 



1. Symmetric modes 



We make use of 
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Eq. (A3), and the symmetries of the square, to obtain from Eqs. (B3) and (B7) the following 
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Subsequently we substitute into Eq. (B13) the decomposition 
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which follows from the Jacobi- Anger expansion e lxcos <t> = J^L_ i n Jn(x)e m ^ and the relation J- n (x) = (—l) n J n (x) 
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Then, using the invariance of the integrand under tp — > it — tp, which implies J^' 2 — > | J*, and performing the 
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substitution cp' = 2ip we obtain 
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The angular integral can be converted into a contour integral over the unit circle in the com plex plane using z = e lip . 
We focus on the case min{n, m} / yfn 2 + m 2 > r) s in which for all 7; the series in Eq. (B15) converges uniformly, so 
that the integral and the sum commute, and the denominator of the first (second) term in the last factor has one real 
root z n< (—z m< ) inside the unit circle and another one z n> (— 5 m >) outside of it given by^l 
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We now evaluate the resulting integral for each term in I by residues to obtain 
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where we have omitted the dependence on qD sin of the 
last factor. Subsequently, we regroup terms in the above 
using Eq. (B17) and substitute Eqs. (B18) and (B16l into 



Eq. (1). Then, we use that here J corresponds to 
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perform the substitution u = c 1 q (for each 7), and inte- 
grate over u). Finally, we eliminate a in favor of cr and 
Pr, introduce the functions fi(x) = ai(x)/(2x 21 ) and 
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and substitute Eq. (Bill, to obtain Eq. (3) specialized 
for n, m odd. 



2. General case 

For the other cases, antisymmetric and mixed- 
symmetry modes, analogous steps as for the symmet- 



ric modes (cf. B12), allow us to obtain from Eqs. (A3 1, 



(B4), (B5), (B8|, and (B9) expressions that differ from 
Eq. ( B13[ ) only in the first two factors of the R.H.S. in- 
tegrand which arc modified via 



cos 2 ( 9 ;)cos 2 ( 9 ;) -> sin 2 (^)sin 2 ( g ;) 
for n, m even (antisymmetric modes) and 
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cos 2 (q' x ) cos 2 (q' y ) -> cos 2 {q' x ) sin 2 (q' y ) 
for n odd and m even (mixed symmetry modes), where 
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q' x = ^sinflcos^ and q' y — sin 6 sin p. The results 



for the other mixed symmetry case [Eqs. ( |B6[ ), (BIO l] can 
be obtained by interchanging the harmonic indices. Sub- 
sequently, we express these factors [Eqs. dB20 1, (B21 1] as 
trigonometric series in p using the Jacobi-Anger expan- 
sion and obtain again Eq. ( |B15[ ), with 
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in the antisymmetric case, and 



a (x) = 1 - J (V2x) , 

ai(x) = —4:J 2 i{x) for /odd, 

ai(x) = -2(-l) l/2 J 2 i(V2x) for I > 0, even 
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in the case of mixed symmetry. Henceforth, the same 



steps followed in Sec. [BTj [cf. Eqs. fB15| , p^ - fB19| ] 
lead to Eq. (3). 

The series in I, corresponding to the last dimension- 
less factor in Eq. (3), can be evaluated numerically using 
Eqs. (|A9])-( |A12| ). Without loss of generality, we assume 
n < m and find that, for ratios m/n not too large, the 
condition J7 7 -\/l + m 2 /n 2 <C 1 implies that prevalently 
77^(1 +m 2 /n 2 )/4 ensuring rapid conver- 
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Equation |l]| is also valid for degenerate resonator modes 
u' R provided that: (i) the degeneracy (possibly split by 
disorderjis associated to a symmetry also satisfied by the 
support 14 or (ii) mode-mixing induced by disorder and/or 
the support is negligible. 

For our structures this gap was < 200nm^C D. 
The average thickness of the membrane is inferred from 
modelling the isotropic etching process in each material. 
The approximate scaling for Qni/Qoi differs from the re- 
sults of applying Eq. |2| by at most a factor of 2 for 
< n < 15 in the relevant regime r\ t < 0.3 and v — 1/3. 
For such modes, ?7 7 £nm «C V n + 1 eventually fails and for 
VjCnm 3> \n 2 — 1/4| one finds that 1/Q„ m exhibits oscilla- 
tory behavior dominated by the SAW contribution. 
In the opposite case min{n, m} / y/ri 2 + m 2 < r] a there are 
one or more types of modes 7 for which at least one of 
the denominators in the last factor in Eq. ( B15 1 has a pair 



of complex conjugate roots on the unit circle, so that this 
factor diverges for certain angles <f> and the convergence of 
the series is no longer uniform, rendering the decomposi- 
tion used invalid. 



